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5 ■ Abstract 

It is shown that the ground state in the Coulomb three-body pdt ion is bound. This ion consists 
of the positevely charged deuterium d and tritum t nuclei and one negatively charged antirpoton 
p. The pdt ion has only one bound S{L = 0)— state which is weakly-bound. The properties 
of this weakly-bound state are investigated with the use of the results of recent highly accurate 
computations. Very likely, the actual proparties of the pdt ion will be different from the results of 
our predictions due to additional contributions from strong interactions between particles. 
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I. INTRODUCTION 



The boundness of the Coulomb three-body system with unit charges formed by the tri- 
ton t (or tritium nucleus), deuteron d (or deuterium nucleus) and one negatively charged 
antiproton p was discussed in a number of our earlier studies jl|, It was shown that this 
system is certainly bound. However, it was also found that the the pdt ion has only one 
bound (ground) S{L = 0)— state and this state is, in fact, a weakly-bound state. Accord- 
ing to the definition of weakly-bound states the ratio of the binding energy e to the total 
energy E of each of these states is very small, i.e. r = -| ^ 1. For weakly-bound states in 
Coulomb few-body systems the dimensionless parameter r must be less than 0.01 (or 1 %). 
A number of Coulomb three-body systems (ions) with unit charges have some weakly-bound 
states. For instance, the P*{L = 1)— states (excited P— states) in the ddfi and dtfi ions are 
extremely weakly-bound. However, the ground S{L = 0)— states are weakly-bound only in 
a very few such systems. In particular, the ground S{L = 0)— state in the three-body pdt 
ion is weakly-bound. It makes the pdt ion unique among all Coulomb three-body systems 
with unit charges, since it contains three heavy particles, but has very small binding energy. 

In general, each weakly-bound state has a number of extraordinary properties. Our 
main goal in this study is to investigate such properties of the ground state of the pdt 
ion. First of all, we want to investigate the properties which are 'typical' for each weakly- 
bound state in the Coulomb three-body systems with unit charges. As follows from the 
results of numerical computations of the weakly-bound states many of such properties are 
substantially mass dependent. Even small variations of particle masses produce noticiable 
changes in the expectation values of large number of properties. To avoid this problem one 
must use extremely accurate variational wave functions. Briefly, this means that the total 
and binding energies of such a state must be determined to the maximal numerical accuracy 
which can be achieved in modern bound state calculations. To achieve this goal we need to 
solve the corresponding Schrodinger equation for bound states in the Coulomb three-body 
system pdt to very high numerical accuracy. The highly accurate wave function \E' obtained 
during this procedure is used later to compute various bound state properties of the pdt ion. 
Some of these properties are of great interest in numerous applications involving the pdt ion. 

It is clear that the structure of the bound state (s) in the pdt ion cannot be explained 
accurately by ignoring contributions from stong interactions between antiproton, deuteron 
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and triton [3|. Indeed, the actual interparticle distances in the ground S{L = 0)— state of 
the pdt ion are only in ^ 50 times larger than the effective radius of the nucleon-nucleon 
interactions (or A^A^— interactions, for short). Therefore, one can expect that the strong 
components of all (three) interparticle potentials in the pdt ion can change the computed 
expectation values, i.e. the bound state properties. In this study, however, we will ignore all 
possible contributions from the strong components of all interparticle potentials. Our main 
goal below is to perform highly accurate analysis of the Coulomb three-body system pdt 
with unit charges. In the next study we are planning to include more realistic interaction 
potentials in our analysis. 



II. THE HAMILTONIAN AND ITS REDUCTION FOR WEAKLY-BOUND SYS- 
TEMS 

In the non-relativistic approximation the Hamiltonian of the three-body pdt ion takes the 
form 4 1 

f^P^2 I rnp^2 I I Qcigte^ ^ QdQpe^ ^ qtqpe^ 

2mp^md mt ^ r^t r^p np 

where Vj = ^) ^) and i = d,t,p. In Eq.([T]) the notation h stands the reduced 

Planck constant, i.e. ^ = and e is the elementary electric charge. For the pdt ion it 

is very convenient to perform all bound state calculations in proton-atomic units in which 

h = l,mp = nip = 1 and e = 1. Here and everywhere below in this study we assume 

that the masses of proton and antiproton exactly equal to each other. The speed of light 

c in the proton-atomic units is c = a"^, where a = |^ is the fine structure constant. In 

proton-atomic units the same Hamiltonian, Eq.([T]), is written in the form 

k-V^ + -V? + v|) + - - -L (2) 



2 ^nid ' ' nit ' ' ' ' Tdt rap np 

where the nuclear masses and mt of the deuterium and tritium nuclei must be expressed 
in terms of the antiproton mass mp which exactly coincides with the proton mass irip. 

The highly accurate wave function of the pdt ion is obtained during numerical so- 
lution of the non-relativistic Schrodinger equation for Coulomb three-body pdt system 
H'^{Yd,Yt,Yp) = E ■ \l'(r(i, rj, Tp), where H is from Eq.(|2]) and < is the total energy 
of this ion. Then by using this highly accurate wave function we can compute a number 



3 



of different expectation values. These expectation values are considered as the bound state 
properties of the ground (bound) state in the pdt ion. 

The Hamiltonian of the weakly-bound Coulomb three-body systems, Eq.(l2]), can be re- 
duced to the sum of the two separated Hamiltonians, i.e. H = Hi + Ho, where the Hamil- 
tonian Hi of the central (or compact) sub-system tp and Ho is the Hamiltonian of the 
outermost particle. For the pdt ion these two Hamiltonians are (in proton-atomic units) 

H. = -l{^V',+Vi)-^ (3) 
2^mt rtp 

and 

rdt rap 

The interaction potential — ^ Hamiltonian Ho, Eq.(jll), is the short-range potential, 

which is, in fact, a polarization potential. This polarization potential acts between the 
positively charged deuterium nucleus (electric charge equals +1) and neutral tp quasi- atom. 
In this two-body picture the deuterium nucleus produces 'sufficiently large' polarization of 
the tp quasi-atom. The interaction between the ion and polarized tp quasi-atom leads 
to the formation of the weakly-bound state in the pdt ion. The boundness of such an ion 
can accurately be evaluated with the use of perturbation theory (see below). Note also 

that the actual, three-particle potential difference between 

the attractive and repulsive Coulomb two-body potentials. All terms in such a potential 
which have non-zero long range asymptotics are compensate each other. Therefore, at large 
distances between the central cluster pt and deuterium nucleus d'^ such a potential almost 
equals to zero. This allows one to reduce the original, very complex, three-body problem to 
the interaction of the central, two-body cluster and one deuterium nucleus, i.e. to the two- 
body problem. Briefly, we can say that the deuterium nucleus produces electric polarization 
of the central cluster. If such a polarization is relatively large, then one (or even a few) 
bound state arise which are stable against dissociation. 

The Schrodinger equation = E-"^ with the Hamiltonian given by Eq.Q is solved ex- 
actly and gives us the hydrogen-like wave functions ip{rtp) of the two-body quasi-atom tp (or 
pt). Then, by using these hydrogenic wave functions we can compute the expectation value 
of the Ho operator, Eq.dl]), which is an operator in respect to the (Cartesian) coordinates of 
the deuterium nucleus d. The origin of the Cartesian coordinates of the deuterium nucleus 
coincides with the center of mass of the tp quasi-atom. The kinetic energy of the deuterium 



nucleus in Eq.(jl]) does not change, if we assume that all wave functions ipivtp) (bound state 
wave functions) have the unit norm. The expectation value of the sum of the two Coulomb 
potentials in Eq.(|l]) can be replaced by some effective potential V{T(i), which in the first 
order approximations is the central potnetial, i.e. V{Y(i) = V{rd). The explicit form of this 
effective potential can be obtained with the use of the perturbation theory. This problem 
is considered in the next Section. It is clear that such a picture is only an approximation, 
since we have ignored a number of additional factors, e.g., contribution from the hydrogenic 
wave functions of continuous spectrum. Nevertheless, the overall accuracy of this model is 
surprisingly high and it still often used for many three-body systems with unit charges. 

III. POLARIZATION POTENTIAL 

As shown above, by using the short-range polarization potential in Eq.(jl]) one can reduce 
analysis of the original three-body Coulomb system with unit charges to the two 'equivalent' 
two-body problems, Eqs.dH]) - dl]). In reality, such an 'equivalent' replacement is very difficult 
to complete, since there are many relations between the explicit form(s) of the model two- 
body potential in Eq.(jl]) and its properties which can be observed in actual, three-body 
experiements. For instnace, as follows from the general theory of the two-body bound states 
in the non-relativistic potential field one finds for the total number N{i) of bound states 
with the angular momentum i 



where the central potential V{r) is the polarization potential mentioned above. For the 
pdt ion we have i = and the mass m in the Bargmann inequality, Eq.([5]), (see, e.g., sl 
and references therein) is the mass of the deuterium nucleus. A slightly more complicated 
analysis indicates that the mass m in Eq.(l5]) must be chosen as the reduced mass of the 
deuterium d and tp quasi-atom. Finally, for the pdt ion one finds from Eq.(l5]) (in proton- 
atomic units) 



Note that all particle masses in the pdt ion are comparable to each other. Therefore, the 
distribution of the electric charge density in this ion is spherically symmetric (or almost 
spherically symmetric). The first order correction upon the charge-dipole interaction be- 




(5) 
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tween the deuterium nucleus d'^ and tp system is zero, since the spatial distribution of 
electric charge in the central tp cluster is spherically symmetric. The second order pertur- 
bation theory leads to the following general form for the interaction V{r) potential between 
the deuterium nucleus and neutral quasi-atom tp in the pdt ion 



The four constants A, a, B and b in this potential must be in agreement with the predic- 
tion which follows from Eq.dH]) and with other experimental data obtained, e.g., from the 
scattering of the deuterium nucleus at the pt system (at different energies). In addition, 
the two-body system with the interaction potential V{r) must have the same (or almost the 
same) bound state properties as the pdt ion. 

It should be mentioned that the two-body polarization potential V{r), Eq.(I7]), carefully 
reconstructed with the use of the known scattering data and bound state properties of the 
pdt ion is of some interest in a large number of applications. For instance, the total number of 
bound states and their approximate geometrical and dynamical properties can be obtained 
with the use of such a potential. On the other hand, such a potential is only approximation 
to the actual three-body potential. Therefore, it is hard to expect that all bound state 
properties of the pdt ion determined with the model two-body potential will be in good 
numerical agreement with the actual properties. 

IV. VARIATIONAL CALCULATIONS 

In general, the highly accurate computations of weakly-bound states in Coulomb three- 
body systems with unit charges are difficult to perform. The main reason for this is obvious, 
since all traditional variational expansions contains 'pieces' which describe contributions 
of fragments from the unbound spectra of two-body systems. In actual computations this 
leads to very slow convergence rate at large dimensions, i.e. when large number(s) of basis 
functions are used in computations. Traditionally, the highly accurate variational com- 
putations of bound states in Coulomb three-body systems with unit charges are performed 
with the use of the exponential variational expansion in perimetric/relative coordinates. The 
explicit form of such an expansion for S{L = 0)-states is 



^ = Ci(j)i{r32, ^31,^21) exp{-aiUi - - 'jiU^) exp{t6iUi + ieiU2 + ifiU^) (8) 




(7) 



N 



1=1 li 
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where Ci are the hnear (or variational) parameters, ai, l3i,'ji,Si,ei and fi are the non-hnear 
parameters and t is the imaginary unit, the function (f)i{r^2, '"si, ^"21) is the polynomial (usu- 
ally quaratic) function of the three relative coordinates Vij. The coefficients of this function 
are fixed and never varied in calculations. The notations Ui,U2 and U3 in Eq.([H]) are the 
three perimetric coordinates: Ui = ^{rij + — Tj^). It can be shown that three periemetic 
coordinates are independent of each other and each of them varies between and +00. 

The variational expansion, Eq.(|8]), provides extremely high accuracy in the bound state 
computations of arbitrary three-body systems (for more discussions, see, e.g., [6|). For highly 
accurate calculations of the ground state in the pdt ion we can assume that all non-linear 
parameters Si,ei,fi equal zero idential, i.e. 5, = 0, Cj = 0,/j = for i = 1,...,A^. This 
means that all varied non-linear parameters in the trial wave function, Eq.(l8]), are real. This 
allows one to re- write the formula, Eq.([8]), to the form 

N 

= ^^Ci exp(-air3i - /J^rsi - 7ir2i) (9) 

i=l £1 

which is called the exponential variational expansion in the relative coordinates r32,r3i and 
r2i. In general, the total energy of the ground state of the pdt ion uniformly depends upon 
the total number of basis functions A^, Eq.Q, used in calculations. 

Now, the phenomenon of slow convergence at large dimensions can be illustrated with 
the use of the asymptotic formula for the E{N) dependence 

m) = E{oo) + ^ + ^ (10) 

at large A^. The formula, Eq. ffTOl) . includes the four parameters E{oo), Ai, A2 and 6. The 
asymptotic energy E{oo) is the improved total energy which is closer to the 'exact' answer 
than the computed E{Ni) values. For bound states in the Coulomb three-body systems 
with unit charges, which are not weakly-bound, the parameter 6 in Eq. ffTOj) does not change 
noticeably with A^. For weakly-bound states in such systems the parameter 6 at large N 
becomes A^-dependent. In many actual cases it decreases when A^ increases. 

To perform numerical calculations in this work we have used another approach devaloped 
in our paper 6^] (see also [7]). In that approach all actual non-linear parameters of the 
method (usually 28 - 40 non-linear parameters) are carefully optimized at some relatively 
large dimension, e.g., for A^ = 1800 and for A^ = 2200. Then the total number of basis 
functions has been increased to A^ = 3500 - 3840. The optimized values of the non- linear 



parameters have not been changed during this (last) step. This simple method allows one 
to obtain results which are significantly more accurate than it is possible to achieve by using 
old-fashion numerical procedures with the same numbers of basis functions. Finally, by 
using our highly accurate wave functions we can determine the expectation values of various 
bound state properties. This problem is discussed in the next Section. 

V. BOUND STATE PROPERTIES 

The total energy and some other bound state properties of the pdt ion can be found in 
Tables I and II. Table I contains the total bound state energies E and the expectation values 
of all three inter-particle delta-functions (5(ri2)), (5(ri3)), (5(r23)) determined in proton- 
atomic units h = l,e = 1 and rrip = 1. Table II includes some other expectation values 
of the bound state properties of this ion. These bound state properties are determined as 
the expectation values of the corresponding operators, e.g. for the operator X we write its 
expectation value X: 



where X — rfj, where n — -2, 1, 2, 3, 4: , X — S{rij),X — 5(ry)5(rjfc), etc. Here and below we 
use the cyclic notations k) — (1, 2, 3) which are very convenient for three-body systems. 
Note that all distances in the pdt ion are in Mp fa 1836 times shorter than the electron- 
nuclear distances in regular hydrogen atom and hydrogen ion. This means that the distance 
between the atniproton and tritium/deuterium nuclei in the ground state of the pdt ion are 

2.5 — 2.6 • 10~^^ cm. Briefly, these distances are only in ?s 20 - 30 times larger than the 
corresponding nucleon-nucleon distances in the light few-nucleon nuclei, e.g., in the ^H, ^He 
and ^He nuclei. This indicates clearly that strong interactions between deuteron, triton and 
antiproton can contribute significantly to the structure and properties of the bound states in 
the pdt ion. In particular, such a contribution will be large for all properties which include 
the interparticle delta-functions, i.e. S{rij) values. 

The quality of the computed delta-functions S{rij) can be checked by comparing the 
computed and predicted values of the corresponding cusp-values Uij which are defined by 
the following relations 




(11) 




(12) 



I S{r,j) I 
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where (ij) = (ji) = (12), (13), (23). For few-body systems interacting by the Coulomb and 
Yukawa-type forces the numerical values of Uij can be predicted and they are not equal zero 
identically. In particular, for the Coulomb three-body system one finds: 



^ij = QiQj : (13) 

rrii + rrij 

where qi and qj are the electric charges and rrii and mj are the masses of the particles. 

By using any two of the three interprticle delta-functions we can define the three-particle 
delta-function 6123 = 5(rij) ■ 5(riA;)- The {6123) expectation value equals to the probability 
to detect all three particles at one spatial 'non-relativistic' point, i.e. inside of the volume 
V = a^ah where a = 1^ is the fine-structure constant, while an = is the Bohr radius. 
This expectation value plays an important role for some systems, e.g., it is used to predict 
the one-photon annihilation in the Ps~ ion. Formally, by using the three-particle delta- 
function 6123 one can try to construct the three-particle cusp operator 1^123, which equals to 
the product of the three-particle delta-function 6123 and a second order (partial) derivative 



d 



However, as follows from the general theory [81] the expectation value of such an 



operator is infinite for an arbitrary Coulomb three-body system. 

The expectaition values of the spatial momenta are defined as follows: 

{^ij) = J J J '^iri2,ri3,r23)r'^'^{ri2,ri3,r23)ri2ri3r23dri2dri3dr23 (14) 

This integral is reduced to the form of the one-dimensional integral of the product of the 
one-particle density matrix and r^. To simplify the following transformations below we shall 
restrict ourselves to the consideration of the bound S{L = 0)— states only. In this case, the 
wave function \l/ is the function of the three scalar interparticle coordinates ri2,ri3 and r23 
only, i.e. \1/ = ^1/(^12, ri3, r23). Therefore, we can define the one-particle density matrix 

and write the following formula for the (r]^) expectation value 

r+oo 

(r-) = / rr2+V(ri2)c^ri2 (16) 

JO 

The idea of this method of 'spatial moments' developed in the middle of 1960's was simple 
and transparent. Indeed, if we can determine a sufficient number of 'spatial moments', then 
we can reconstruct, in principle, all three one-particle density matrixes p(ri2), p(ri3) and 
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p{r2z)- In turn, this will lead us to useful conclusions about spatial distributions of the two 
and three particles in the system. However, this idea has never been applied to reconstruct 
the actual density matrixes. Furthermore, since the middle of 1960's we can approximate 
the wave function of any three-body system to very high numerical accuracy. By using these 
wave functions we can determine all spatial momenta to extremely high accuracy (see Table 
II). However, a few questions about spatial momenta are still remain unaswered for three- 
body systems. The most interesting question is related to the 'singular momenta', i.e. to the 
(r^) expectation values with m = —3, —4, —5, . . .. These expectation values are singular, 
but their principal and regular parts are needed in some problems, e.g., the values {r^"^) are 
used to evaluate the lowest-order QED correction (or Araki-Zucher correction) for Coulomb 
two-electron atoms and ions (see, e.g., [q] and references therein). Our Table II contains the 
epectation values of the regular part of the (r^^^) expectation value which is designated as 
{r~j^)fi. These expectation values are related to the corresponding principal parts ((r^^^)) by 



the following relations (see, e.g., {q] and references therein) (r^j^) = {r~f')R + 47r((5(rjj)). 

Table II also contains expectation values which are used to measure various interparticle 
correlations in three-body systems. The spatial interparticle correlations are evaluated with 
the use of the three cosme-functions r: 

= C-^^) (17) 

where a • b means the scalar product of the two vectors (a and b), r^j = rjj, but Vij = —rji, 
since = rj — r^. The sum of the three Tj^ values always exceeds unity, i.e. it can be 
written in the form 

Ti2 + ri3 + r23 = 1 + 4 ■ / (18) 

where the /—value is 

/ = I I *)=/// I *(r„.n3. r,3) r ^i^n.n,rMudr., 



2J J J \ '^{ui,U2,U3) l"^ UiU2U3dUidU2du3 (19) 



In some earlier works the expectation values of the three scalar products Tij ■ were 
proposed to be used to describe spatial interparticle correlations in the three-body systems. 
Note that the three interparticle vectors rjj,rjfc,rfcj always form a triangle. Therefore, we 
can write 

Tki + Yij = Ykj , or ■ Tij = ^{rl + rfj - rlj) (20) 
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In other words, the expectation values of all three scalar products Tij ■ r^fc can be expressed 
from the known (rfj) values. This means that these scalar products are not independent 
properties. Analogous situation can be found for the three scalar products pj • of each 
two single particle momenta. It follows from the conservarion of the total momentum P 
and its three components (P)j. Indeed, in the center-of-mass system we can write P = 
Pi + P2 + Ps = 0. From here one finds three following identities 

P.-Pj = ^(p'+P?-Pfc) (21) 

The same identities can be written for the corresponding expectation values. Therefore, 
the three expectation values of the scalar products pj ■ pj are uniformly related to the 
single particle 'kinetic energies' \{pf) and cannot be used to describe dynamical correlations 
between three point particles as independent values. 



VI. ANNIHILATION AND FUSION RATES 



In reality, the three-body pdt ion is not stable, since it decays either by the proton- 
antiproton annihilation (at rest), or by the nuclear (d, t)— fusion. Note that the proton- 
antiproton annihilation in the pdt ion proceeds as the tritium-antiproton annihilation, or 
deuterium-antiproton annihilation. The corresponding equations are: 

pdt = 7i^ + n+ + n- + d + 2n + Qi{MeV) (22) 

pdt = 7r^ + Kt + KL + t + n + Q2{MeV) (23) 

and 

pdt = p + '^Re + n+ n.lMeV (24) 

respectively. Here the notation n stands for the neutron, while vr^ and designate the 
corresponding tt— mesons (charged and/or neutral), and denotes the i^— mesons, 
while "^He means the alpha-particle. In general, the proton- ant iproton (or triton-antiproton 
and deuteron-antiproton) annihilation can proceed with the use of other reaction channels 

n 

(see, e.g., and references therein). The energy released during antiproton annihilation 



is usually around hundreds and even thousands of MeV [10|] (this depends upon the an- 
nihilation channel). In the fusion reaction, Eq.f l24p . the antiproton p plays the role of a 
particle-'catalyzator' which simply brings the two positive particles closer to each other. 
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The probability of formation of the bound, two-body p^He system (also called pa-system) 
after such a reaction in the pdt ion in is relatively large (~ 5 % according to our approximate 
evaluations) . 

Currently, the both tritium-antiproton and deuterium-antiproton annihilations at small 
energies of the colliding particles are not well studied phenomena. In particular, there 
are substantial deviations between experimental data and theoretical predictions based on 
different theoretical models. Moreover, the role of antiproton annihilation at neutron(s) 
has never been evaluated properly. By using our expectation values of the interparticle 
delta-functions determined for the ground state of the pdt ion (see Table I) we can evaluate 
the relative contributions of the antiproton annihilation at the triton and deuteron nuclei, 
respectively. Indeed, for the antiproton annihilation rates (F) at the deuteron and/or triton 
we can write the following formulas: 

rpd = {cipp + Ctpn) i^pd) , rp = (cLpp + '^Ctpn) {^pt) (25) 

where the notations p, p, n stand for the proton, antiproton and neutron, respectively. From 
here one finds 

Fp ^ + (M ~ 7.8161771198656- f|il) (26) 

where ^ = ^ is the ratio of the corresponding 'basic' annihilation rates (the rates of the 
neutron-antiproton and proton- antiproton annihilation at zero interparticle distances). It is 
extremely interesting to evaluate the ratio ^ experimentally. This will give us an improtant 
information about the antiproton-neutron annihilation at rest. Also, it follows from Table I 
for the pdt ion the antiproton annihilation rate at the triton is ~ 3.9084 times more likely 
than annihilation at the deuteron. 



VII. HYPERFINE STRUCTURE 



The hyperfine structure of the pdt ion is of great experimental interest, since the magnetic 
moments of the proton and antiproton are equal to each other by their absolute values, but 
have opposite sings. Note also that all three particles in the pdt ion have non-zero magnetic 
moments (or spins). Therefore, we can expect that the hyperfine structure of the ground 
state of the pdt ion will be sufficiently reach. Moreover, such a structure and corresponding 
hyperfine structure splittings can easily be observed in modern experiments. However, the 
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determined values of the hyperfine splittings of the pdt allows us to understand a large 
number of interesting details related to the long-range asymptotics of the strong interparticle 
interactions. 

For the bound S{L = 0)— state in the pdt ion the hyperfine structure is determined by 
solving the eigenvalue problem for the Hamiltonian {AH)h.s. which is responsible for the 
spin-spin (or hyperfine) interactions. The general formula for such a hyperfine Hamiltonian 
{AH)h.s. in the case of the pdt ion is written as the sum of the three following terms. Each 
of these terms is proportional to the product of the factor and the expectation value 
of the corresponding (interparticle) delta-function. The third (additional) factor contains 
the corresponding factors (or gyromagnetic ratios) and scalar product of the two spin 
vectors. For instance, for the pdfi ion this formula takes the form (in atomic units h = l,e = 
l,me = 1) (see, e.g., I^], 11|) 

/Am 27r 29iP9d,^, x , 27r 29p9t,., ^w ^ 



^2vr^,^ (27) 



„2 



where a = 1^ is the fine structure constant and rap = nip is the proton mass. Note also that 
the presence of the common factor ^ in all three terms in Eq. (|2711 indicates clearly that the 
proton-atomic units mentioned above are more appropriate for this problem. Nevertheless, 
below we shall continue to use the usual atomic units ^ = 1, e = 1 and rrie = 1. To recalculate 
the expectation values of the 5— functions from proton-atomic to atomic units one needs to 
use the factor ^ 6.1905094020072 ■ 10^ The expression for {AH)h,s., Eq.(I27]), is, in fact, 
an operator in the total spin space which has the dimension {2sp + l)(2sd + ^){'2st + l) = 12. 
In our calculations we have used the following numerical values for the constants and factors 
in Eq.(l27D: a = 7.297352586 ■ 10"^ and rrip = 1836.1527012me. The ^-factors for the 
antiproton, deuteron and triton are deteremined from the formulas: g„ = ^^^-.Qd = 4^ and 
Qt = where Aia and are the corresponding magnetic moments and spins of the particle 
a. For the triton, deuteron and antiproton we have /j = |,Jd = 1 and Ip = ^, while the 
magnetic moments of these particles (in nuclear magnetons) are Mp = —2.792847386, Aid = 
0.857438230 and Mt = 2.97896247745. 

The hyperfine structure and hyperfine structure splittings of the ground state of the pdt 
ion can be found in Table III. Traditionally, these values are expressed in Mega Hertz, or 



13 



MHz. To-recalculate our results from atomic units to MHz we used the conversion factor 
6.57968392061 -10^ MHz/a.u. 12|. As follows from Table III the hyperfine structure of 
the pdt ion include twelve levels which are separated into the four following groups: (1) the 
group of five spin states with J = 2, (2) the upper group of three states with J = 1, (3) one 
state with J = and (4) the lower group of three states with J = 1. The energy differences 
between these four groups of states are: 1.112524269-10^ MHz, 6.853155051-109 MHz and 
5.977402022-10'' MHz. These values can be measured quite accurately by using modern 
experimental technique developed to measure atomic hyperfine splittings. 

Note again that such a hyperfine structure corresponds to the three-body pdt ion with the 
pure Coulomb interaction potentials between particles. In reality, corrections to the strong 
interactions between particles will change the expectation values of all interparticle delta- 
functions. Therefore, we can expect some changes in the hyperfine structure and hyperfine 
structure splittings. Such corrections to the hypefine structure splitting of the pdt ion will 
be the goal of our next study. 



VIII. CONCLUSION 



We have determined the total energy and some bound state properties of the weakly- 
bound (ground) S{L = 0)— state in the pdt ion. It is shown that this weakly-bound state has 
a 'two-body' cluster structure, which is represented as the motion of the deuterium nucleus 
d~^ in the central field of the 'central' pt cluster. The interaction between the deuterium 
nucleus d'^ and the two-body pt cluster is the regular dipole-charge interaction, which is 
often called and considered as the polarization potential. The use of the two-body cluster 
model and approximate reconstruction of the polarization potential allows one to evaluate 
qualitatively many bound state properties of the pdt ion. In highly accurate calculations 
of the ground state of the pdt ion performed for this study we have used our advanced 
optimization procedure developed in ^ and later improved in [7]. This strategy works 
perfectly for all Coulomb three-body systems, including weakly-bound and cluster systems. 

A number of bound state properties of the pdt ion have been determined to very high nu- 
merical accuracy with the use of our bound state wave functions. In particular, by using the 
expectation values of the delta-functions we have investigated the hyperfine structure of the 
ground S{L = 0)— state in the pdt ion. In this ion there are twelve levels of hyperfine struc- 
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ture which are separated into four different groups. The energy differences between these 
groups of states are A12 = 1.112524269-10^ MHz, A23 = 6.853155051-109 MHz and A34 = 
5.977402022-10'' MHz. These values are called the hyperfine structure splittings. The ob- 
tained values of the hyperfine structure splittings must be compared with the corresponding 
experimental values. 

Another interesting question for this ion is related to a direct comparison of the rates of 
the (d, t)— nuclear fusion and rates of the antiproton annihilation at the deuteron and tri- 
ton, respectively. These problems have solutions which are very sensitive to the contribution 
from strong interparticle interactions. Therefore, the comparison of the future experimental 
results and theoretical predictions will lead us to a large number of improtant conclusions 
about long-range asymptotics of the strong interactions. Currently, we have developed a 
number of effective numerical methods for such computations. Moreover, as follows from 
the results of fisrt calculations performed with the use of 'realistic' interparticle potentials 
(Coulomb potential plus term which describes strong interactions) we can say that the con- 
tributions from the antiproton-triton and antiproton-deuteron strong interactions play the 
leading role. Analogous contribution from the deuteron-triton realistic potential is signifi- 
cantly smaller. Therefore, to make accurate predictions for the ground state of the pdt ion 
we need to know the interaction potential between the antiproton and triton and deuteron, 
respectively. The contribution from the triton-deuteron strong interaction is also improtant, 
but it is not crucial. In general, the accurate reconstruction of the potential of strong 
interactions require a detailed knowledge of all its scalar, spin-orbital (L - S) and tensor 
components. To reconstruct such a potential for each pair of interacting particles one needs 
to know a large amount of scattering data |l^ and additional information about all bound 
two-body states. 
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TABLE I: The convergence of the total energies E and expectation values of the two-particle delta- 
functions for the ground (bound) S{L = 0)— state of the pdt ion (in proton-atomic units). A'^ is the 
total number used in calculations. 

N E {631) {Sptl (Spd) {Sd^ 

3300 -0.38119 08996 43549 79944 1.035725171596-10-^ 2.65020902768-10-2 1.97194857-10-'^ 
3500 -0.38119 08996 43549 80319 1.035725171777-10"^ 2.65020902581-10"^ 1.97194865-10"^ 
3700 -0.38119 08996 43549 80574 1.035725171405-10"^ 2.65020903462-10"^ 1.97194853-10"^ 
3840 -0.38119 08996 43549 80725 1.035725171481-10"^ 2.65020902513-10"^ 1.97194857-10"^ 



17 



TABLE II: Some bound state properties of the ground (bound) S{L = 0)— state of the pdt ion (in 
proton-atomic units). The particle 1 means deuteron, 2 desinates triton, while 3 stands for the 
antiproton. 



{X) 


(ij) = (32) 


(ij) = (31) 


(ij) = (21) 




0.6452584980114 


0.3265829532244 


0.20945965194864 


{{nkrjkY 


^) 0.090672197685358 0.13833503600526 


0.17327231991926 


ink") 


0.894094641294119 0.289502249538164 0.0635359440628808 


{nj) 


2.535959239261 


5.822464028824 


6.5332990137423 


(4) 


9.369595134746 


52.35601818235 


57.621835007682 


(4) 


46.65710362569 


645.7495736644 


676.62342378591 


(4) 


294.4990484287 


10212.70615551 


10290.298340574 




0.76584205188631 


0.34062934917079 


0.0766784190714 




0.084729933577141 


0.05594048174487 


1.217745967-10^^ 




-0.7496066938104 


-0.66655644341827 


1.1986352684990 


(6) 

v-- 


-0.7496066901228 


-0.66655635263091 


1.1986366724385 




0.069022844172409 0.22440893908996 


0.27170237010133 


{pi ■ Pi) 


0.006716068397607 0.01060650083845 


0.012940402539131 


(rife • Tj-fe) 


50.30412902764317 7.31770598003839 


2.05188915470772 




-0.4103977567344 


-1.8132102882933 


0.02747944473455 



(")The values of /, ((ri2ri3r23) ^) and ((J123), respectively. 
(''^The exact values of the two-body cusps. 
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TABLE III: The levels of hyperfine structure e and hyperfine structure splittings A of the ground 
bound S{L = 0)— state of the pdt ion (in MHz). 



ej A 

pdt{J = 2) 2.3642257710-109 

pdt{J = 1) 2.2529733441-109 1.112524269-10^ 
pdt{J = 0) -4.6001817067-10^ 6.853155051-10^ 
pdt{J = 1) -4.6599557269-10^ 5.977402022-10^ 
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